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The stability of a metal thin films on a dielectric substrate is conditioned by the magnitude of 
the interactive forces at the interface. In the case of a non-reactive interface and weak adhesion, 
the minimization of free surface energy gives rise to an instability of the thin film. In order to 
study these effects, Pt thin films with a thickness of 50 nm were deposited via ion-beam sputtering 
on yttria stabilized zirconia single crystals. All Pt films were subjected to heat treatments up to 
973 K for 2 h. The morphological evolution of Pt thin films has been investigated by means of 
scanning electron microscopy (SEM), atomic force microscopy (AFM) and standard image analysis 
techniques. Three main observations have been made: i) the deposition method has a direct impact 
on the morphological evolution of the film during annealing. Instead of hole formation, that is 
typically observed as response to a thermal treatment, anisotropic pyramidal shaped hillocks are 
formed on top of the film, ii) It is shown by comparing the hillocks' aspect ratio with finite element 
method (FEM) simulations that the hillock formation can be assigned to a stress relaxation process 
inside the thin film, iii) By measuring the equilibrium shapes and the shape fluctuations of the 
formed Pt hillocks the anisotropy of the step free energy and its stiffness have been derived in 
addition to the anisotropic kink energy of the hillock's edges. 

PACS numbers: 68.60.Dv, 68.55.J-, 81.16.Rf 
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I. INTRODUCTION 

Thin metal films on dielectric substrates are thermo- 
dynamically instable. Their stability has been subjected 
to research under several aspects: thermodynamics and 
kineticsi - — , mass transport via surface diffusioni^, 
impact of surface energy anisotropics 6 9 , fingering 
instabilities ^ 1Q i 11 , Ostwald ripening of islands^ 2 -, hole 
pattern 13 and hillock formation^"— .Most of the funda- 
mental theoretical work has been carried out by Srolovitz 
and Safran who developed a complete stability theory for 
thin films covering kinetics 20 and energetics 2 ^. 
In the case of strained layers, the shape instability 
leads either to equilibrium-shaped hole or to equilibrium 
shaped hillock formation^ 2 -"— depending on the com- 
peting relaxation mechanism. While the formation of 
hillocks as a consequence of stress relaxation is observed 
for various thin film materials 1 ^ 19 ! 23 ! 25 " — , the formation 
of equilibrium shaped holes has been affirmed recently 
both experimentally and by kinetic Monte Carlo simula- 
tions 6 . However, the nature and transition between these 
two competing instability mechanisms is far from being 
fully understood. This is mostly due to a lack of experi- 
mental data. 

The main objective of this paper is to establish a rela- 
tion between the observable macroscopical changes dur- 
ing hillock formation on a strained thin film and its un- 
derlying configurational forces, e.g., the kink energy. The 
investigation will focus on how the presence of an internal 
stress field in the thin film triggers the thin film insta- 
bility caused by hillock formation. In a second step, it 
is addressed, how the equilibrium hillock shapes can be 



used to determine critical stability-related quantities like 
the step line stiffness /? and the kink energy e£. There- 
fore standard finite element modeling (FEM) has been 
chosen alongside the analysis of the anisotropic hillock 
shape fluctuations, which is successfully applied to de- 
termine the step energies of 2D islands during thin film 
growth 2 ^"— . 

In the present study, Pt thin films on single crystalline yt- 
tria stabilized zirconia (Zr0 2 ) have been chosen as model 
metal/ceramic systems in order to study the hillock for- 
mation. In equivalence to Au/Zr0 2 32 , the Pt thin film 
can be regarded as strained and its interface to Zr0 2 
as semicoherent 33 , due to a lattice parameter misfit of 
e m = (af - a s )/a s = 0.31. 

This paper is structured as follows: Section [TT] derives the 
basic theoretical models necessary to describe the forma- 
tion of hillocks. Section IIIII deals with the experimental 
framework. In Section HV( detailed results are presented 
and discussed. The final section [V] encompasses a sum- 
mary of the findings and conclusions. 



II. THEORETICAL BACKGROUND 
A. Basics 

During deposition of metal thin films via sputtering, 
the kinetic energy E^ n of the deposited atoms generally 
exceeds their thermal energy i^therm- For the deposited 
film, this results in a metastable configuration that tends 
to equilibrate, once subjected to temperature by anneal- 
ing, Joule's heating or radiation. The morphological sta- 
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FIG. 1. Illustration of the experimentally observed regimes of 
hillock formation (a) during thermal treatment of a dense flat 
film (b) hillocks form on the film surface with a regular hexag- 
onal shape (c) due to further annealing holes and secondary 
hillocks form in the vicinty of the primary hillock 



bility of a metal thin film on a dielectric material is thus 
conditioned by the aspect ratio, the interaction across 
the interface 34 and the tendency of the thin film to re- 
duce its free energy, e.g., due to stress relaxation. The 
way in which the reduction of the free energy manifests 
in the evolution of the film morphology depends strongly 
on the competing relaxation mechanisms. While in typi- 
cal thin film agglomeration scenarios, defect related local 
perturbations cause a film rupture and a decrease of sur- 
face areai, the formation of hillocks on the contrary is 
the direct response to a delocalized strain field in the 
thin film. The origin of this strain field is attributed to 
a lattice mismatch e m or a growth stress induced by the 
deposition technique 35 . It is noteworthy, that in contrast 
to thin film agglomeration the formation of the hillock is 
usually accompanied by the increase of the surface area 
as shown schematically in Fig. HJ The change of the total 
free energy AF of a uniformly strained film with volume 
V to a (partially) relaxed film with hillock on top can be 
expressed in terms of the change in strain energy density 
AW and surface energy A<£, 

AF = AW + A$ = (W - Wo) + (S - * ). (1) 

Whereby A<£ > 0, if only hillocks are formed. Hence, 
the reduction of the free energy can only be caused by 
a decrease in elastic energy which scales with the initial 
elastic energy 



W = VMe 2 n 



(2) 



where M 



27T/1 



with the Poisson ratio v and the shear 



modulus [i. The relaxed strain energy W throughout the 
film volume V is considered to be equal to the induced 



change in the elastic strain field 
hillocks. Thus 



by the formation of 
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CijklCijCkldV, 



(3) 



where Cijki is the stiffness tensor of the material. In gen- 
eral, Eq. [3] has to be solved numerically. 

In order to facilitate the calculation, the three dimen- 
sional hillock can be approximated by a two dimensional 
frustum with rotational symmetry, as shown in Figure [2j 
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FIG. 2. Schematic drawing of an elastically strained film in- 
cluding a hillock in the shape of a frustum, rotationally sym- 
metric around the axis M, with the top radius Ri, the bottom 
radius R2 and the height h 



By doing so, the calculation of the change in surface en- 
ergy 



A^> = jfQ — 7/^0 



(4) 



is simplified as it is characterized only by the uniform 
surface energy 7/ and the surface area ft given by the 
three geometrical parameters Ri, R2 and h. The surface 
after a hillock has formed reads, 



n = 7r(R 1 +R 2 )-jRi - 2R X R 2 



Rl 



■h 2 +7rRl + (n -7rR 2 2 ). 

(5) 

W obviously depends on the film morphology and there- 
fore on Ri,R 2 and h. fto is the surface area of the initially 
flat film. 



B. Hillock Shape Analysis 

Hillock formation as well as thin film agglomeration re- 
quire the motion and creation of atomic steps. Thereby 
the general anisotropic line stiffness /3 of the atomic steps, 
serves as the key parameter to investigate and under- 
stand the morphological evolution of a wide class of thin 
film instabilities 36 . In analogy to islands growt h 29 ! 37 , the 
threefold symmetric hillocks can be treated as two dimen- 
sional shapes and their contour line can be expressed in 
terms of a Fourier-series with the general form 



Rm(0) =R + J2 a i- sin(n,(# - 0*)). 



(6) 



i=l 



Whereby Ro, a$, 6i serve as fitting parameters. The 
step free energy per unit length /?, or step line tension, 
is related to the hillock shape R^t by the 2D Wulff con- 
struction 24 , which has been proven by Burton et al 38 to 
be 



m = x 



R it 



fit 



(7) 



Thereby <p = 6 — arctan (deRnt/Rfu) denotes the normal 
to the equilibrium shape for each R^ t (0). It is notewor- 
thy, that Eq. [71 establishes the proportionality between 
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the measurable 2D shape Rfa(0) and its free energy f3. 
For small shape fluctuations, the elongation of the step 
contour line has to be taken into account and the step 
line tension is replaced by the step line stiffness = 
f3((j)) + d^/?. Similar to Eq. [7J the step line stiffness (3 is 
orientation-depended and related to the hillock's curva- 
ture k by 



0(61) = A. 



(R 2 



■fit 



R lt 
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A 
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(8) 



Once the equilibrium shape of a hillock is reached, the 
curvature n of steps along the densely packed directions 
is zero and /? = /?. Conversely, all shape fluctuations 
cause a curvature k of the step. Using the Terrace Step 
Kink (TSK) model 31 , this curvature is related to the ther- 
mally activated formation of kinks along the step with the 
length of n atomic units an. In the case of an unrestricted 

TSK-model, the step line stiffness j3 is given by 



P = (2 an kb T / a\) sinh 2 



2fa,T' 



(9) 



where €k is the required kink formation energy and 
an = 0.277 nm and a\ = 0.240 nm are the unit lattice 
spacing parallel and orthogonal to the step edge, respec- 
tively. 

The only variable in the upper equations, which is not 
yet determined, is A. Kodambaka et al£^ established a 
generalized formulation for anisotropic 2D crystal shapes, 
which directly relates A and the experimentally accessible 
hillock shape fluctuation function g. The hillock shape 
fluctuations g are determined by calculating the devia- 
tion of the experimentally measured hillock contour from 
the fitted equilibrium shape given by Eq.[6]and is defined 
as follows, 



g(0) = (r exp - Rm)/RfLt- 



(10) 



In order to relate g to A, two fluctuation sensitive func- 
tions 



x(0) 



y/Rj t + 2d 9 Rl t - ifetdeeifct 



and 



p(0) =g-Rm 



(11) 



(12) 



are defined. The analytical functions of x an d p are ex- 
pressed in terms of a Fourier series with Fourier coeffi- 
cients Xn and p n . By applying the equipartion theorem, 
which states that in equilibrium every DOF has the same 
energy (E) 



\hbT, it can be shown that 



Kg) 



Xn\ 



(13) 



where A/" max = 0.5 • N atoms ^Q is defined as one half of 
all atoms along the contour line dQ of the hillock. A 
detailed reproduction is given in reference 39 . 



III. EXPERIMENTAL 



A. Sample Preparation and Characterization 



The Pt/Zr0 2 (single-crystalline) system is immiscible 
and characterized by a lattice mismatch of e m = (aj — 
a s )/a s = 0.31. A chemical inert interface is formed for 
all temperatures below 1273 K— . Pt layers of 50 nm 
in thickness were deposited at room temperature by ion 
beam sputtering (pbase = 1 • 10 -7 mbar) upon the single- 
crystalline substrates that were pre-cleaned using iso- 
propanol. Before deposition, the single crystalline sub- 
strates were cleaned in the Ar-ion beam of the sputter- 
ing chamber for 5 s. The substrates coated with the 
thin platinum film were annealed in a muffel furnace 
for 2 h at 923 K and 973 K. The annealing tempera- 
tures are well below the melting temperature of platinum 
Tm — 2042 K, hence volume diffusion of Pt in Pt is pro- 
hibited. The morphology of the samples was studied via 
high resolution AFM, using a Mobile S (Nanosurf), and 
SEM (Zeiss Leo 1530). For all acquired hillocks, the as- 
pect ratio a and the equilibrium shape r{0) have been 
determined. 



B. Finite Element Modelling 



In order to confirm the role of incompatible strains due 
to a lattice mismatch as the primary reason for hillock 
formation, a thermoelastic FEM model using COMSOL 
Multiphysics has been developed. The film is assumed 
to be subjected to an extensional mismatch strain e m . 
The origin of this strain is the large lattice mismatch / 
between film and substrate which can be treated as a 
thermoelastic deformation. A nonuniform elastic strain 
€ij is created by assigning a thermal expansion strain 
equal to — e m in the thin film. A frustum is chosen to 
represent the hillock shape in first approximation, see 
Fig.[2j The frustum is rotationally symmetric around the 
axis M and features therefore three-dimensional effects. 
In analogy to the AFM- analysis, the frustum's bottom 
radius R2 and the hillock height h define the hillock's 
aspect ratio a = I1/R2. The film is treated to be isotropic, 
hence the initial strain energy is Mf e m V where Mf is the 
biaxial modulus of the film and V the total film volume. 
The total free energy change AF of the system during 
hillock formation is simulated starting from a flat film 
with increasing hillock height h and coevally decreasing 
bottom radius R2. The material-specific constants used 
in the FEM model are listed in Tab. U 



TABLE I. FEM model parameters for the elastic constants 
E, v and the surface energies 7pt, 7ysz and 7pt/ysz 



Parameter 


Value 


Lit. 


E 


1.63 • 10 11 Pa 


41 


V 


0.390 


41 


7Pt,m 


1.656 


42 


7YSZ 


1.927 


43 


7Pt/YSZ 


1.2 


1 



IV. RESULTS AND DISCUSSION 
A. Basics 
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(a) 923 K (b) 973 K 



FIG. 3. SEM images of different stages of hillock formation 
on Pt(lll) with corresponding annealing temperature, (a) 
Hillock formation on top of the dense film, (b) Additional 
hillock growth and holes forming in the vicinity of the hillock. 



For all isochronically annealed (t = 2 h) Pt thin films, 
the stages of hillock formation have been analyzed. In 
Fig. [3j the two different stages of hillock formation on 
the Pt thin film with an initial thickness h = 50 nm 
are shown representatively. For sufficiently low anneal- 
ing temperatures T a , the film relaxation is solely driven 
by the formation of hillocks while all other regions on 
the film are smooth, see Fig. |3(a) . Energetically speak- 



ing the gain in free energy is dominated by the strain 
relaxation, hence W < The observed hillock shapes 
in Fig. 3(a) are highly anisotropic and characterized by 
a truncated triangle with a threefold symmetry. This 
symmetry is defined by the different edge lengths (A, B) 
of the hillock which entail the energetical differences of 
the series of steps forming these edges. The base of the 
hillock is enlarged in respect to its height h, which re- 
flects the expected kinetics for hillock growth described 
by Tersoff and Tromp22. In contrast to recent molecu- 
lar dynamics simulations of hillock growth^, the hillock 
nucleation in this work does not demand a seed grain po- 
sitioned on top of the film. 

By increasing T a , the morphological evolution of the thin 
film undergoes a transition indicated by the coexistience 
of two different relaxation modes: the additional hillock 
formation and film rupture in the vicinity of the hillock 
[Fig. |3(b)] . Thereby the film rupture along the long edges 
(B-steps) of the hillock corresponds to the onset of sur- 
face roughening in regions of local stress enhancement, as 
shown in Fig. |3(b)| However, the additional aggregation 
of hillocks preferentially occurs at the short edges of the 
hillock, which coincides with the low energy steps (A- 
steps). While the hillocks base is significantly enlarged, 
the mean hillock height is not severely affected by the 
heat treatment. These observations are in good agree- 
ment with the predictions for shape changes induced by 
strain relaxation made by Tersoff 22 . 
The hillock shapes and their aspect ratio a have been ac- 
quired by high resolution atomic force microscopy, with 
a lateral resolution < 0.2 nm. In total 35 hillocks were 
analyzed. In Fig. 4(a) a typical 80 nm high hillock is 
shown. The hillock has a characteristic threefold sym- 
metry, the edges of the hillock are steep. The base ra- 



dius R2 and top radius R\ of the hillock differ slightly. 
For all hillocks, formed after annealing at T a = 923 K, 
the contour line of the base has been measured and the 
center O has been calculated. Furthermore the radii ta 
and vb correspondi ng to the A and B steps have been 
determined, see Fig |4(b)| Due to the different orienta- 
tion of the steps A( 10 o) and £?(m), the step formation 
energies are generally different. The anisotropy mani- 
fests in the ratio of step formation energies, which results 
from the Wulff-relation /3 a/ ' Pb = ta/vb- The Wulff re- 
lation follows from Eq. provided that dgR = 0, and 
thus applies only to the extrema of the hillock radius r. 
The mean ratio of the step free energy resulting from the 
measured A and B radii is (3 a/ (3 b = 1.56(11) and given 
in Tab. [TTJ The found anisotropy and energetically pre- 




FIG. 4. (a) 3D AFM image (750 • 750nm 2 ) of a hillock 
on Pt(lll) annealed at 923 K. (b) Hillock perimeter of the 
formed hillock obtained via standard image analysis tech- 
niques, the differing radii r a and rb indicate a anisotropy in 
the line tensions for A- and B-steps. 

ferred formation of 5-step is in accordance with both 
growth experiments on pit formation on Pt(lll) single 
crystals 37 and ah initio calculations on Pt(lll) step en- 
ergies 44 , see also Tab. [Ill However, the difference in the 
measured step energy ratio in this work and the ones of 
Michely et al 37 and Boisvertj 4 ^ is significant and can be 
related to an unequal mechanism that causes the relax- 
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(a) a = 0.19 



(b) a = 0.39 



(c) a = 0.70 



FIG. 5. (a-c) Strain energy density balance calculated using a thermoelastic FEM- model of a 50 nm Pt thin film on Zr0 2 for 
rotationally symmetric hillocks with increasing aspect ratio a. 



TABLE II. Comparison of the determined length ratio ta/tb 
of the hillock, and the step energy ratio Pa/Pb with data 
from literature 



Parameter 


Value 


Lit. 


Pa/Pb 


1.56(11) 


this work 


Pa/Pb 


1.15(2) 


exp. 37 


Pa/Pb 


1.13 


theor. 44 



the facilitated formation of dislocations and mass trans- 
port caused by the enhanced stress in the film 23 . 
The relaxation process has been modeled by FEM on 
the grounds of Eqs. [2][3] for aspect ratios a ranging from 
0.1 to 2 for four different lattice mismatches e m ranging 
from the theoretical lattice mismatch e m = 0.31 to a mis- 



ation process. While in the present work the relaxation is 
controlled by a stress field generated from the lattice mis- 
match between film and substrate, the formation of pits 
with anisotropic shape found in literatur e) 3744 originates 
from a 2D to 3D growth transition, which is facilitated 
by a large tensile stress of the Pt(lll)-surface. 
This detailed impact on hillock formation of a mismatch 
induced stress field is addressed in the next section. 



B. Finite Element Modelling 




-O- s m =0.075 [33] 

A e m =0.10 
-O-e m =0.15 
— n— e =0.31 [theor. ej 



0.5 1.0 
aspect ratio a 
(a) 



1.5 



In consequence of the previous findings we assume 
that growth of hillocks is caused by a relaxation pro- 
cess whereby the minimization of the free energy AF 
is dominated by strain relaxation. Using the concepts 
derived in Sec. Ill Al the experimentally revealed aspect 
ratio a = h/R\ obtained by AFM are compared using 
FEM modeling with the predicted optimal trade-off be- 
tween hillock shape and newly formed surface area Q. 
The total film volume is conserved, thus the measured 
effect is solely defined by a redistribution of mass. In 
Fig. [5j the strain energy balance Wo — W for three differ- 
ent aspect ratios a is shown. In comparison to its initial 
value Wo, the strain energy density W is decreased sig- 
nificantly (blue regions) by the formed hillock for all a 
in Figure Furthermore with increasing a an increasing 
strain energy density in the vicinity of the hillock base 
has been observed. As shown in Fig. 2J these regions of 
high strain energy density coincide with regions of film 
rupture after further thermal treatment. This is due to 




simulated (FEM, Sm =0.075) 
aspect ratio (AFM, hillocks) 



0.3 0.6 0.9 

aspect ratio a 

(b) 



FIG. 6. (a) Total free energy change as function of the aspect 
ratio a as determined from the FEM-model. (b) Evolution of 
the total free energy as function of the aspect ratio a (FEM) 
including the experimentally determined aspect ratio distri- 
bution of hillocks (AFM). 
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FIG. 7. (a) AFM image of a hexagonal shaped hillock on Pt(lll) film annealed at 923 K (b) Measured equilibrium hillock 
shape plotted in polar coordinates radius R vs. angle (open circles) fitted with Eq. [6] (solid line) (c) Polar plots of the fitted 
equilibrium shape R(0) and the step free energy /3(<t>) with A = 1 



match e m = 0.075 resulting from HRTEM analysis^ 3 -. It 
has to be noted, that the implemented strain in the film 
does not necessarily result from a lattice mismatch e m but 
can also be due to an intrinsic growth strain. The total 
free en ergy change AF as function of a and e m is shown 
in Fig. 6(a)| AF(a) resembles a Lennard- Jones potential, 
whose potential- well depth scales with e m . However, the 
position of the minimal free energy stays unaffected from 
changes in e m at a = 0.20, hence the found relation can 
be regarded as general scaling law under the condition 
that \AW\ > |A$|. 

In order to validate the functional relationship between 
AF and a, the experimentally quantified distrib ution of 
hillock aspect ratios a has been compared in Fig. |6(b)| to 
the minimal free energy configuration predicted by the 
FEM modeling. Despite the geometrical simplification 
in the FEM model, a good agreement between the mea- 
sured mean hillock aspect ratio a = 0.19(2) and the pre- 
dicted minimal free energy configuration at a = 0.20 of 
the thin film has been found. Thus, it can be concluded 
that the found shape change of the thin film is domi- 
nated by a strain relaxation process, whereby the formed 
hillocks correspond to the minimal energy configuration 
of the film. 



C. Hillock Shape Analysis 

In order to provide a deeper insight into the fun- 
damentals of hillock formation, the shape fluctuation 
method 29 ^"— is used to determine the kink energy e 
and the corresponding line stiffness /?. The anisotropic 
kink energy e is a physical key property that impacts 
the physics and decay kinetics of two-dimensional Ost- 
wald ripening 29 , island growt h 37 ! 39 , vacancy island for- 
mation 45 , fingering instabilities 7 -, thin film agglomera- 
tion^ and hillock formation. 

Due to the anisotropy of the found equilibrium hillock 



shapes, e and j3 depend on the step orientation A, B. By 
using the unrestricted TSK model (Eq. [9j), the kink en- 
ergy can be determined from the line stiffness /?, which is 
directly related to the chemical potential of the hillock 31 . 
In the following discussion, the presented results apply 
exclusively to hillocks formed at T = 923 K. At this 
temperature a hillock is not restricted or influenced by 
another hillock or hole in its vicinity. The shape fluctu- 
ations of hillocks with radii ranging from 150 to 550 nm 
have been analyzed. For each hillock, the equilibrium 
shape R and g were determined from the AFM data. A 
typical data set is shown in Fig. It is co mpos ed of 
the initially measured hillock by AFM (Fig. |7(a)[ ), the 
educe d pol ar-plot of the hillock perimeter fitted by Eq. [6] 
(Fig. |7(b)[ ) and a combined mapping of the calculated 
step free energy (Eq.[7]) with A = 1 and the hillock shape 
in Cartesian coordinates (Fig. |7(c) ). Due to entropy in- 
duced shape fluctuations 30 , there is a measurable differ- 
ence between the fit-function and experimental data in 
Fig. |7(b)| As long as the shape fluctuations g are small, 
typically < 10%, the fluctuations of the step free energy 
5j3 are small and hence can be approximated by the step 
line stiffness f3 (see Eq. [9]). This condition is complied 
with all hillocks, as the chosen fit-function (Eq. [6]) agrees 
well with the experimental hillock shapes measured by 
AFM (R 2 > 0.92). 

The shape fluctuation function g{6) has been determined 
for each hillock from the AFM data. Subsequently g(6) 
and Rftt have been inserted into Eqs. [11] and [12] to cal- 
culate and p(0). The two functions have been ex- 
panded into Fourier series. The resulting Fourier coeffi- 
cients Xn and p n are used to obtain A for each hillock. 
Thereon the step line stiffness (3 has been calculated nu- 
merically for the A and B steps using Eq. [8] The kink 
energies ca,b result from solving Eq. [9] 
The average values of the step line stiffness /3a,b an d the 
kink energies €a,b are given and compared to literature 
data in Tab. IIII1 The experimental step line stiffness of 



7 



TABLE III. Comparison of the experimentally determined 
stiffness clu/Ba and an /3b and the kink energies ca and cb all 
in meV with data from literature 





a\\f3 A 


cl\\Pb 


€A 


€B 


(this work) 


335(50) 


1582(220) 


186(30) 


312(60) 


Exp. 48 


482(32) 


1471(102) 


144(3) 


206(4) 


Exp. 49 






167 


167 


Num.^ 






180 


250 



A and B steps obtained at T = 923 K shows a significant 
anisotropy and is in the magnitude of stiffnesses obtained 
from vacancy island formed on Pt(lll) single crystals at 
T = 713 K 45 . Thereby /3a is slightly smaller compared to 
the equivalent step line stiffness measured by Ikonomov 
et al 45 and /3b slightly larger. This difference obviously 
refers to the enhanced anisotropy of the A and B step 
length in this work and to the different annealing tem- 
perature. It is assumed, that the found hillock shape 
represents the equilibrium shape in respect to the domi- 
nant configurational forces in this thin film geometry at 
T = 923 K. Although the experimental framework in 
this work differs significantly from the one of Ikonomov 
et a/4 5 -, a remarkable agreement of the step line stiffness 
is achieved. 

The same applies to the calculated kink energies ca,b, 
which are in good agreement with ab-initio calculations 50 
as well as with the experimental results of Ikonomov 45 , 
see Tab. IIIII In contrast to the data of Giesen et a/. 49 , 
the kink energy e is significantly step- type dependent. 
The numerical values of the measured kink energies cor- 
respond well to the theoretical results of Feibelman 5 ^ at 
T = OK and therefore cast the postulated temperature 
dependence of the kink energy 51 into doubt. 
However, although the absolute values for ca,b found in 
this work are slightly larger than in reference 45 , the ratios 
of the kink energies are with S € = ^ « 0.5 in good accor- 
dance. The numerical values of ca,b can now be used to 
determine the step depended bond energy Ja,b = 2ca,b 
which serves as a key parameter in solid-on-solid KMC 
models 7 simulating thin film instabilities. It becomes ob- 
vious that if one wants to foresee and understand the 
fundamentals of thin film instabilities, it is inevitable to 
consider the general anisotropic character of the kink en- 
ergy which directly affects the motion of atomic steps. 



minimization of the strain density AW and the other to 
the minimization of the surface energy A<£. These two 
competing contributions to the total free energy cause a 
transition from pure hillock formation at T < 923 K to a 
coexistence of hillock formation and film rupture in the 
vicinity of the hillock edges at T > 973 K. The observed 
competition of relaxation mechanisms in this work is in 
excellent agreement with the predictions made by Ter- 

These findings are substantiated by the performed 
FEM simulations that clearly indicate that formation of 
hillocks on strained films cause a minimization of the 
thin films free energy as function of the hillocks aspect 
ratio a. The predicted aspect ratio that corresponds to 
the minimal free energy configuration of the thin film 
is in excellent agreement with experimentally measured 
mean aspect ratio. Therefore it is assumed that the ob- 
served hillocks possess equilibrium shape. In addition, it 
has been shown by FEM that the local maximum of the 
stress field in the thin film coincides with experimentally 
observed regions of rupture. 

The induced hillock formation is governed by the motion 
of atomic steps on the surface. The analysis of the hillock 
shapes by high resolution AFM revealed quantitatively 
two key parameters of the motion of atomic steps: the 
anisotropic step line stiffness (3a,b and the anisotropic 
kink energy ea,b- The found values are in good agree- 
ment with both experimental 45 and theoretical 5 - - data 
from literature. On closer inspection, it can be seen that 
the hillock shape anisotropy found in this work slightly 
differs from the observed anisotropics in formed vacancy 
islands on Pt(lll)^ 5 - or grown Pt-islands 52 . This is due 
to the different state of the surface and the presence of a 
substrate material in case of a thin film geometry. 
It is noteworthy, that the choice of the deposition tech- 
nique, ion beam sputtering in this work and magnetron 
sputtering in a previous work 1 , severely impacts the 
pathway of the thermal instability of the thin film. While 
in a previous study thermal annealing of the thin film 
led to hole formation at defect associated perturbations, 
the identical thermal treatment caused the formation of 
hillocks in the present work. Hence the inherent thin film 
properties like the distribution of grain boundary ener- 
gies or the internal stress induced by the deposition tech- 
nique have a decisive impact on the competing instability 
mechanisms and their subsequent shape changes. 



V. CONCLUSIONS 

In essence, it has been shown that the formation of 
hillocks in ion beam sputtered Pt thin films is a result of 
a complex stress relaxation process. The kinetics of the 
inherent mechanisms exhibit a different thermal activa- 
tion and depend differently on the stress field inside the 
thin film. It is shown that two in general independent 
physical processes control the morphological evolution 
and kinetics of hillock formation: one is attributed the to 
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